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Abstract 

We study the long-time behavior of the probabihty density Qt of the first exit time from a 
bounded interval [—L, L] for a stochastic non-Markovian process h[t) describing fluctuations at a 
given point of a two-dimensional, infinite in both directions Gaussian interface. We show that Qt 
decays when t ^ oo as a power-law where a is non-universal and proportional to the ratio 

of the thermal energy and the elastic energy of a fluctuation of size L. The fact that a appears to 
be dependent on L, which is rather unusual, implies that the number of existing moments of Qt 
depends on the size of the window [—L, L\. A moment of an arbitrary order n, as a function of L, 
exists for sufficiently small L, diverges when L approaches a certain threshold value L„, and does 
not exist for L > Ln- For L > Li, the probability density Qt is normalizable but does not have 
moments. 

PACS numbers: 02.50.-r, 05.40.-a, 05.70.Np 
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Non-equilibrium surface growth and dynamics of interfaces belong to an important mul- 
tidisciplinary branch of non-equilibrium statistical mechanics, which received much interest 
within several last decades due to its numerous practical applications. To name just a few, 
we mention such diverse fields as crystal growth, molecular beam epitaxy, fluctuating steps 
of metal surfaces, evolution of liquid-liquid or liquid- vapor interfaces, and growing bacterial 



colonies 
equations 




. A number of discrete atomistic models 



|2|, y, |j] and stochastic evolution 
has been proposed, revealing a generic scale invariance as manifest 
in the power law behavior of the interface width and correlation functions of the interface 
height. Much effort has been also expended in understanding universal aspects of different 
theoretical models and naturally occuring processes 12, 13||. 

While earlier works have focused mainly on scaling behavior, more recent effort has been 
concentrated on evaluation of the distribution functions of characteristic properties of equi- 
librium or non-equilibrium interfaces. For different models of interfaces several probability 
distribution functions have been determined, including, e.g., the distribution of the width of 

it distribution in one HI] and two- dimensions 



leig hts in the steady state [10], maximal heig^ 



12 1 , distribution of spatially averaged height 131] or of the height at any fixed point in space 
141 ]. and the density of local maxima or minima of heights 15|. Last but not least, statis- 



tics of persistence, i.e., the probability that the interface height at a given point remains 
persistently above (or below) its initial value during some time interval, has been studied 



3 



17|. 



In this paper we address a basic extreme value problem for a non-Markovian stochastic 
process h{t), where h{t) is a height (measured relative to the averaged value) at a fixed 
point of a non-equilibrium two-dimensional, infinite in both directions Gaussian interface, 
separating two phases in ci = 3. Using the approach of Refs. l3], we determine the long-time 
asymptotical behavior of the "survival" probability Pt that h{t) does not leave within the 
time interval (0,t) a bounded interval [—L,L], and correspondingly, define the tail of the 
probability density Qt of the first-exit-time t from [—L,L], Qt ~ —dPt/dt. We show that 



18 



19|, and 



in a striking contrast to the one- dimensional case, where Pt ~ exp 
hence, Qt has all moments, for a two-dimensional Gaussian interface the probability Qt is 
characterized by a power-law tail of the form 



Q{t) 



1 



t ^ K ^ exp 



(1) 
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T being the temperature (measured in the units of the Boltzmann constant ks) and k, - the 
interfacial tension. The exponent a in Eq.([T]) is non-universah 

where C is a constant, such that 7r/16 < C < 7t/8. The exponent a is thus proportional 
to the ratio of the thermal energy and the elastic energy of a fluctuation of size L. As 
a consequence of the L-dependence of a, not all moments of Qt exist and, what is rather 
unusual, the very number of existing moments of Qt depends on the size of the window 
[-L,L]. 

Consider a two-dimensional interface whose local heights hn^m (t), measured relative to 



the averaged value, obey an infinite set of Langevin equations [20|: 
where — oo < n,m < oo and Ct"'™^ are independent Gaussian white-noise processes: 



In Eqs.© and (jl]), ^ is the friction coefficient, the bar denotes averaging over thermal 
histories, 6n,n' is the Kronecker-delta symbol and 6{t) - the delta-function. In what follows 
we set, for simplicity, ^ = 1, such that the appropriate dimensionless "time" variable will be 
K,t. Dependence on ^ can be restored in our final results by a mere replacement T T/^ 
and K, k/C,- 

Note that Eqs.(l3]) describe, in particular, the time evolution of a spatially discretized 
Edwards- Wilkinson interface 5], as well as model A Langevin dynamics [2^ of the Weeks 
columnar model {2]]] or of a coarse-grained interface in a three-dimensional Ising model above 



the roughening temperature 22 1. 



We suppose that initially the interface is fiat and hn^m{t = 0) = for all n and m. 
Applying to Eqs.Q a discrete Fourier transformation, solving the resulting equation and 
inverting the solution, we find that, for a given thermal history, the local height at the 
origin h{t) = hofi{t) is defined as a portfolio of an infinite number of weighted independent 
Gaussian processes: 

/t 00 00 

dre''-- E E Ini2KT)Im{2Kr)Ct'r\ (5) 
n=— 00 m=— 00 



where /„(2Kr) is the modified Bessel function of order n. From Eq.([5]) one finds that at 



sufficiently long times 



T 



h?it)^- — \n{Kt). (6) 
^ ' Ann ^ 

Our aim is now to define, for a non-Markovian process h{t) in Eq.Q, the probability Pt 
that h{t) has not ever crossed the boundaries of the interval [—L, L] within the time interval 
{0,t) given the initial condition h{t = 0) = 0, i.e., Pt = Prob(max|/i(t)| < L\h{0) = 0), where 
ma.x\h{t)\ is the largest absolute value achieved by h{t) within the time interval (0,t). Once 
Pt is determined, we will get an access to the behavior of another important probability - the 
first-exit-time probability density Qt = Prob(t' > t\h{0) = 0), where t' = min{r|/i(r) = ±L} 
is the time when h{t) first hits either of the boundaries; hence, Qt is defined as 

dP 

Qtdt = --^dt = Prob {t<t'<t + dt\h{0) = 0) . (7) 

Focusing on the large-t asymptotical behavior, we note that it is not really important 
how we define Ci"'*"^ - as continuous in time functions or as discrete processes, provided that 
we keep all essential features of noise. We thus divide, at fixed t, the interval (0,t) into 
small subintervals A, (such that AA^ = t), and assume that 

/■(n,m) _ f^^\ ci^'^) /'o^ 

where [x] denotes the fioor function. In other words, we assume that (^^"''"•' is constant and 
equal to a/2T/A S'^"'™'' within the k-th subinterval, k = 0,1, . . . , N — 1, where {5'^"'"^''} is an 
infinite set of independent identically distributed random variables with normal distribution 
N[0,1]. 

Then, h{t) can be expressed as a weighted sum of an infinite number of independent 
discrete noise processes: 

oo oo N 
?i=— oo n=— oo 1=1 

with weights 

/ rp \ 1/2 p2KAl 

a;"-)= -— / due-'-lMUu). (10) 

Squaring Eq.Q and averaging the resulting expression with respect to the distribution of 
i.i.d. variables {Sj^'"^^}, we get 

TV 



h^{N) = J2^f, (11) 



1=1 



which expression introduces an effective variance af, defined as 



^' = E E G^^'O (12) 

^ " ' ^ ' ^'^<^ '2kA(/-1) J2kA(/-1) 



2 2" 



n=— oo m=— oo 

where 

Noticing that exp(— 2a;)/Q (x) is a monotonically decreasing function of x, one finds that af 
is bounded by two monotonically decreasing functions of /: 

~ 2 

exp (-8kA/) II (4kA/) < < exp (-8/tA(/ - 1)) (4/tA(/ - 1)) . (14) 

When / ^ 1, these bounds become very sharp and hence, with a good accuracy, 

d] ^ 2AT exp (-SkA/) ll (4kA/) . (15) 

Substituting Eq. (fT5|l into Eg. (ITT!) , we recover the result in Eq.(l6]). 

Define now the following event: An A^-step discrete-time trajectory h{N), Eq.Q, com- 
mencing at the origin, does not leave the interval [—L,L], or in other words, the maximal 
absolute value, max|/i(A^)|, which process h{N) achieves for a given realization of noise, is 
less than L. Probability of such an event we denote as P/v = Prob(max|/i(A^)| < L). 

Now, condition max\h{N)\ < L implies that the absolute value of any ascending partial 
sum 

oo CO N 

hdN) = E E E ^^'"''stiK (16) 

n=—oo m=— oo /=A''— A;+l 

which define the values of the local height h{N) at consecutive discrete "time" moments k, 
k = 1,2, . . . , N, is bounded from above by L. 

We prefer, however, to deal with the descending partial sums: 

oo oo fc 

K{N)= E T.^t'"'^SN-i\ (17) 

n=— oo m=—OD 1=1 

which define the trajectory hk{N) evolving in the inverse time N — k and shifted by a 
constant (realization-dependent) value hN{N). Clearly, P/v = Prob(max|/i(A^)| < L) = 
Prob(max|/i'(A^)| < L). 

Let I(max|/i'(i(:)| < L) be the indicator function: 

{1 , max\h'(N)\ < L, 
(18) 
0, max\h'{N)\ > L. 



and Rl{x) - a rectangular function, such that: 



Rl(x) 



dy si\i{Ly) 



,iyx 



1, 


\x\ < L, 


< 1/2, 


X = ±L, 


.0' 


\x\ > L. 



(19) 



Then, using the definition of the descending partial sums in Eq. ([T71) and Eq. ([T^ . we write 
down Eg. (1181) as the following iV-fold integral: 



l(max\h'{N)\ < = JJ l(^|/i'^(A^)| < L 

k=l 



n 

k=i 



dyk sm{Ly, 



IT 



Vk 



k)^iykK{N)_ ^20) 



Averaging the indicator function in Eg. (1201) with respect to the distribution of i.i.d. random 
variables {5'^"''^''}, and changing integration variables (see Refs.|l8| for more details), we 
find eventually the following exact representation of P/v: 



N 



L N 



dhk 



exp 



2^! 



(21) 



where is fixed, = 0. Once P/v is known, the desired survival probability Pt can be 
obtained as an appropriate limit: Pt = limA-*o,N ~,oo Pn , with AA^ = t kept fixed. 

y and one has to 



18|, we construct 



The A^-fold integral in Eg. (12 II) can not be, of course, evaluated exacts 
resort to controllable approximations. Here, using the approach of Refs. 
rigorous lower and upper bounds on P/v in Eg. fl2T|) . taking an advantage of the following 
property of P^ in Eg. ([21]) ui- 

P/v = P/v((Ti, 5-2, (T3, . . . , 5-jv) is a monotonically decreasing function of any variable dk- 
This signifies that replacing any or all 5"^ by S(fc), such that cxk < S(A;), we will decrease 
the right-hand-side of Eg. (l2T]) and arrive at the lower bound on P/v; if, on contrary, we will 
replace one or all dk by S(/c), such that dk > S (/;:), we will increase the right-hand-side of 
Eg. (pT!) and obtain an upper bound on P/v. 

We start with a lower bound on P/v. As we have already noticed, exp(— 2x)/g (x) in 
Eg. (1121) is a monotonically decreasing function of x. Conseguently, setting in the integrand 
^2 = 0, and integrating over du2, we have 

rp r-2K,Ak 
^ J2KA{k-l) 



due~^''ll{u) 



Aupp) _ Aupp) 
^k ^k-l ' 



(22) 
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where 



t 



(upp) 



T 



K 



due~^''ll{u) 







Therefore, P/v in Eg. (1211) is bounded from helow by 

-L 



Pn> 



exp 



{hk - h 



k-l 



(23) 



(24) 



where is fixed, /iq = 0. 

One notices now that the expression on the right-hand-side of Eq. 



probabihty that an A^-step trajectory of a Brownian motion evolving in "time" tu does 



describes the 

{upp) 



not leave an interval [— L, L], which yields, in the leading order. 



Pn > exp 



71" Aupp) 



8L2 



exp 



2kAN 



du exp(— 2'u)/q (m) 



(25) 



Defining the asymptotical large- behavior of the integral in the exponential in Eq. (l25|l and 
turning to the limit A ^ 0, — > oo with AA^ = t kept fixed, we arrive at the following 
lower bound on Pf. 

]-n( PA TT T 

(26) 



ln(P+) n T / ^9 ,r^. 

^ ' > —, Kt > exp(KLVT). 



ln(/tt) - 8 2/tL2 

Consider next an upper bound on P/v, Eq. fl2T]) . Using the inequality in Eq. (fT^ and an 
evident inequality: exp(— 2x)Jq (x) > 1/ {2ttx + 1), which holds for any x > 0, we have 

2AT 



di > 2 AT exp {-8k Ak) (AnAk) > 



SixKAk + 1 



t 



(low) Alow) 



k-l ' 



whith 



t 



(low) 



T 



ln(87rKA(A;+l) + l). 



Attk 

Hence, in virtue of Eq. (|2711 . P/v in Eq. (l2T]) is bounded from above by 



P/v < 



L N 

n 



dhk 



exp 



^=1 J27r(tf'")-t£^) 



{hk — hk-iY 



2 t 



.{low) Alow) 



k-l 



(27) 



(28) 



(29) 



where ho is fixed, /iq = 0. 

One recognizes now that the rhs of Eq. (l29|) is the probability that an A^-step trajectory 
of a Brownian motion evolving in "time" t^'""'-' does not leave an interval [— L, L]. Hence, 
the following upper bound on P/v is valid: 



Pn < exp 



T^Alow) 



exp 



(30) 



which yields the following upper bound on Pt: 

m}^-T62^' ««»-P(-^VT). (31) 

Since Pt is, evidently, a monotonically decreasing function of time, we infer that the large-t 
asymptotical behavior of Pt is described by a power-law of the form Pt ~ (/tt)"", where 
a = CT/2k,L'^ and C is a constant, such that tt/16 < C < tt/S. Consequently, the long-time 
tail of the first-exit-time probability density Qt is a power-law, Qt ~ Eq.([l]). 

In conclusion, we have shown that the probability density Qt of the first-exit-time t from a 
bounded interval [—L, L] for a non-Markovian process h{t) describing fluctuations at a given 
point of a two-dimensional, infinite in both directions Gaussian interface, is characterized by 
a power-law tail in Eq.([T]). The exponent a is non-universal and proportional to the ratio of 
the thermal energy and the elastic energy of a fluctuation of size L, and thus depends on L. 
We note that a power-law behavior of Pt is not counterintuitive, and could, in principle, be 



expected from a heuristic estimate Pt ~ exp(— [23] combined with a logarithmic 
growth of the second moment of h{t), Eq.([6]). On the other hand, the dependence of a on 
L is a rather unusual fact which entails unusual behavior of the moments of Qf, namely, 
the number of exisiting moments of Qt appears to be defined by the size of the window in 
which the stochastic process h{t) is observed. More specifically, the number n of existing 
moments depends on the relation between L and a discrete set of characteristic lengths 
Ln = (CT/2Kn)^/^; the condition that Qt has exactly n moments is, in fact, equivalent to 
the requirement that L obeys the following double-sided inequality: -L„+i < L < L^. In 
other words, a moment of an arbitrary order n, as a function of L, exists for sufficiently 
small L, diverges when L approaches L„, and does not exist for L > Ln- For L > Li, the 
probability density Qt is normalizable but does not have moments. We note finally that an 
analogous behavior can be expected for the interfaces defined by noisy MuUins equation [g] 
in d = 5 or by general linear Langevin equations described in Refs. in d = 3. 
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